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This is intend to provide an overview of the theory and phenomenology parts of the TMD 
(Transverse Momentum Dependent parton distribution and fragmentation functions) 
studies. By comparing with the theoretical framework that we have for the inclusive 
deep inelastic lepton-nucleon scattering and the one-dimensional imaging of the nucleon, 

I try to outline what we need to do in order to construct a comprehensive theoretical 
framework for semi-inclusive reactions and the three dimensional imaging of the nucleon. 

After that, I try to give an overview of what we have already achieved and make an 
outlook for the future. 
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1. Introduction 

With the deep going of the study of the nucleon structure, three dimensional imaging 
has become the frontier and a hot topic in recent years. It is commonly recognized 
that the three dimensional imaging contains much more abundant physics on the 
nucleon structure and the properties of QCD. The study was initially triggered 
by the experimental finding of striking single-spin asymmetries (SSA) in inclusive 
hadron production in hadron-hadron collisions. Gradually it grows into a field aim¬ 
ing at a comprehensive picture of nucleon structure including spin and transverse 
momentum dependences. 

The one dimensional imaging of the nucleon is provided by the parton distribu¬ 
tion functions such as the number densities, q{x), the helicity distributions, Aq(x), 
and the transversities, Sq{x), for quarks of different flavors. In the three dimensional 
case, i.e. where the transverse momentum is also considered, not only the direct ex¬ 
tensions of them to include transverse momenta are involved, but also many other 
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correlation functions such as the Sivers function, the Boer-Mulders function, the 
pretzelocity etc exist. Moreover, higher twist effects become also important and 
need to be considered consistently. 

The study on the three dimensional imaging of the nucleon is in a rapid develop¬ 
ing phase and it is not so easy to present a comprehensive overview of all different 
aspects of studies. Here, 1 choose to do the job in the following way: First, 1 will 
try to make a brief review of what we did in one dimensional case with inclusive 
DIS. In this way, 1 hope that 1 can show you the main line of what we need to do in 
three dimensional case. Then 1 will try to summarize what we have already achieved 
along this line and what we need to do next. For the sake of space, 1 will mainly 
concentrate on the discussions but keep as less equations as possible. An extended 
version is prepared and will be published in a special issue of Frontier of Physics. 

2. Inclusive DIS 8z the One Dimensional Imagining of the Nucleon 

Our studies on the structure of a fast moving nucleon started with the inclusive 
deep inelastic scattering (DIS) e.g. e“ -I- iV e~ X. We recall that, under one 
photon exchange, the differential cross section is given by the Lorentz contraction 
of the known leptonic tensor and the hadronic tensor W^i, {q, p, S) (where p and S 
are the 4-momentum and polarization vector of the nucleon). Information of the 
structure of the nucleon is contained in the hadronic tensor W^i/(q,p, S). 

The theoretical framework for inclusive DIS has been constructed in the following 
steps. First, we studied the kinematics and obtained the general form of the hadronic 
tensor by applying the basic constraints from the general symmetry requirements 
such as Lorentz covariance, gauge invariance, parity conservation and Hermiticity. 
We found out that the hadronic tensor is determined by four independent structure 
functions Fi, F 2 , gi and 52 , where the first two describe the unpolarized case and 
the latter two are needed for polarized cases. 

Our knowledge of one dimensional imaging of the nucleon starts with the “in¬ 
tuitive parton model” that is very nicely formulated e.g. in [I]. Here, it was argued 
that, in a fast moving frame, because of time dilation, quantum fluctuation such as 
vacuum polarizations, can exist quite long. In the infinite momentum frame, such 
fluctuations exist for ever. In this case, a fast moving nucleon can be viewed as a 
beam of free “partons”. The probability of the scattering of an electron with the 
nucleon is equal to that of the scattering with a parton convoluted with the number 
density of the parton in the nucleon, i.e., 

|AI(eiV->■ eA)p = ^ f dxfq{x)\M{eq ^ eq)\^, (1) 

<? 

where fq (x) is the number density of q in the nucleon. In this way, we obtained the 
famous results, F 2 (a:,( 52 ) = 2xFi{x,Q‘^) 

and so on. Here, I would like to point out that, with this intuitive parton model, 
we are doing nothing else but the impulse approximation that we often used in 
describing a collision process where we do the following approximations. 
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• during the interaction of the electron with the parton, interactions between 
the partons are neglected; 

• the scatterings of the electron with different partons are added incoherently; 

• the electron interacts only with one single parton. 

Although the intuitive model is elegant and practical, we are not satisfied since 
it is not easy to control the accuracy. The proper quantum field theoretical (QFT) 
formulation is given by starting with Feynman digram Fig. 1(a) where we obtain, 

1FW(9,p,5) = (2) 

where k is the 4-momentum of the parton, (q, k) = ^+ ^)^v{2T:)5+{{k — qY) 
is a calculable hard part and the matrix element 

= J d‘‘ze*'"^(p,S'|V5(0)V’(z)|p,S'), (3) 

is known as the quark-quark correlator describing the structure of the nucleon. By 
taking the collinear approximation, i.e. taking k « xp, and neglecting the power 
suppressed contributions, we obtain exactly the same result as that obtained using 
Eq. ([1]) based on the intuitive parton model. [Here, we use light-cone coordinate 
kh = (fc+, fcj_) and take h = (1,0, 0 _ l ), n = (0,1,0j_), n± = (0, 0, nj_). Also, we 
choose the coordinate system such that p = p+h.] At the same time, we obtain also 
a QFT operator expression of fq{x), 

fg{^) = I ^e^^P^^-{p\m^Y{z)\p), ( 4 ) 

which is indeed the number density of parton in the nucleon. However, from this 
expression, we see immediately that it is not (local) gauge invariant! To get the gauge 
invariant formulation, we need to take into account the multiple gluon scattering 
given by the diagram series shown in Fig.l(a-c). The contribution from each diagram 
is expressed as a trace of a hard part and a matrix element. E.g., for Fig. 1(b), 

^ j ^^■^^TT[Hl];^Yki,k2,q)(i>l^Hki,k2,p,S)], ( 5 ) 

4'^p\ki,k2,p,S) = J d‘^zd‘^ye'^'^^^~^‘^^^~^^'>^{p,S\'ilj{0)Ap{y)'ipiz)\p,S), ( 6 ) 

where L in the superscript denotes left cut. The matrix element is now a quark-j- 
gluon(s)-quark correlator, where j is the number of gluons. We also see that none 
of such quark-_)-gluon(s)-quark correlators is gauge invariant. 



(a) (b) (c) 


Fig. 1. Examples of the Feynman diagram series with multiple gluon scattering considered for 
'y* + N ^ q-\- X with (a) j = 0, (b) j = 1 and (c) j = 2 gluons exchanged. 
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To get the gauge invariant form, we need to apply the collinear expansion as 
proposed in Refs. [21141 . which is carried out in the following four steps. 

(1) Make Taylor expansions of all hard parts at ki = Xip, e.g., 

{k,q)= TLW {x) + V + • • • , (7) 

where ujp^ is a projection operator defined by uipP = Pp^ — np-n^ . 

(2) Decompose the gluon field into Ap{y) = A'^{y)fip + ojpP Ap'{y). 

(3) Apply the Ward identities such as, 


dH^^Xx) 

dkp 





HjSjxi) 

X 2 — xi — ie 


( 8 ) 


(4) Add all terms with the same hard part together and we obtain Wpi, {q, p, S) = 
J2j c (j = 0; 1; 2, • • • and c is the different cut) where e.g. for j = 0, 


WW [q, p, 5) = T I ^Tr [Hj°J {x) (k, p, ,5)], (9) 

^^°\k,p,S) = J 5 '|^^( 0 )£( 0 ; 2 /)i/'(?/)b, 5 '), ( 10 ) 


where 13(0; y) = l3'l'(oo; 0)13(oo; j/), and £{oo;y) = ^ (( "W) glands 

for path integral), is the well-known gauge link obtained in the collinear expansion. 

In this way, we have constructed the theoretical framework for calculating the 
contributions at the leading order (LO) in pQCD but leading as well as higher twist 
contributions in a systematical way. The results are given in terms of the gauge 
invariant parton distribution and correlation functions (PDFs). We also see that 
the PDFs involved here are all scale independent. This is because we have till now 
considered only the LO pQCD contributions, i.e. the tree diagrams. 

Because the hard parts in wid) ’s such as that given by Eq. (jH]) are only functions 
of the longitudinal component x but independent of other components of the parton 
momentum fc, we can simplify them. E.g., for j = 0, it reduces to. 


WW(g,p, S) = ^lp+dxTr[Hjd}{x) ^)], (H) 

¥°\x,p,s) = J ^e^^p^y-{p,s\mmy-my-)\p,s). (12) 


We see clearly that only one dimensional imaging of the nucleon is relevant in 
inclusive DIS. We also see that the PDFs are defined in terms of QFT operators 
via the quark-quark correlator *5') by expending it in terms of y-matrices 

and the corresponding basic Lorentz covariants. 

To go to higher order of pQCD, we take the loop diagrams that describe gluon 
radiations and so on into account. After proper handling of these contributions, we 
obtain the factorized forirP where the PDFs acquire the scale Q-dependence govern 
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by QCD evolution equations. In practice, PDFs are parameterized and are given in 
the PDF library (PDFlib). 

In summary, for studying one dimensional imaging of the nucleon with inclusive 
DIS, we took the following steps, 

• General symmetry analysis leads to the general form of the hadronic tensor 
in terms of four independent structure functions. 

• Parton model without QCD interaction leads to LO in pQCD and lead¬ 
ing twist results in terms of Q-independent PDFs without (local) gauge 
invariance. 

• Parton model with QCD multiple gluon scattering after collinear expansion 
leads to LO in pQCD, leading and higher twist contributions in terms of 
Q-independent but gauge invariant PDFs. 

• Parton model with QCD multiple gluon scattering and “loop diagram con¬ 
tributions” after collinear approximation, regularization and renormaliza¬ 
tion leads to leading and higher order pQCD, leading twist contributions 
in factorized forms as functions of Q-evolved and gauge invariant PDFs. 

In the following, I will follow this four steps and summarize what we achieved 
in the three dimensional case. Before that, I would like to recall the following two 
of historical developments that may be helpful for to us to construct the theoretical 
framework for the TMD case. 

First, as mentioned, the study of three dimensional imaging of the nucleon was 
triggered by the experimental observation of SSA in the inclusive hadron-hadron 
collision. It was known that pQCD leads to negligibly small asymmetry for the hard 
part but the observed asymmetry can be as large as 40%. The hunting for such large 
asymmetries last for decades with following milestones: 

• In 1991, Sivers introduced^ the asymmetric quark distribution in a trans¬ 
versely polarized nucleon that is now known as the Sivers function. 

• In 1993, Boros, Liang and Meng proposed^ a phenomenological model that 
provides an intuitive physical picture showing that the asymmetry arises 
from the orbital angular momenta of quarks and what they called “surface 
effect” caused by the initial or final state interactions. 

• In 1993, Collins published^ his proof that Sivers function has to vanish. 

• In 2002, Brodsky, Hwang and Schmidt calculated^ SSA for SIDIS using an 
explicit example where they took the orbital angular momentum of quark 
and the multiple gl uon scattering into account. 

• Collins pointed oulP^ that the multiple gluon scattering is contained in 
the gauge link and that the conclusion of his proof in 1993 was incor rect 
because he forgot the gauge link; Belitsky, Ji and Yuan resolved2Il2! the 
problem of defining the gauge link for a TMD parton density in light-cone 
gauge where the gauge potential does not vanish asymptotically. 

Another historical development concerns the azim utha l asymmetry study in 
SIDIS. It was shown by Georgi and Politzer in 1977 thaifl^ final state gluon radia- 
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lions lead to azimuthal asymmetries and could be used as a “clean test to pQCD”. 
However, soon after, in 1978, it was shown by Cahn thalP^ similar asymmetries can 
also be obtained if one includes intrinsic transverse momenta of partons. The latter, 
now named as Cahn effect, though power suppressed i.e. higher twist, can be quite 
significant and can not be neglected since the values of the asymmetries themselves 
are usually not very large. 

The lessens that we learned from these histories are in particular the following 
two points, i.e., when studying TMDs, 

• it is important to take the gauge link into account; 

• higher twist effects can be important. 

Both of them demand that, to describe SIDIS in terms of TMDs, we need the 
proper QFT formulation rather than the intuitive parton model. 


3. TMDs Defined via Quark-Quark Correlator 

The TMD PDFs of quarks are defined via the quark-quark correlator given 
by Eq. (HOI) (after integration over k~). A systematical study has been given in 
Ref. [IS] and a very comprehensive treatment can also be found in Ref.[T2]. Here, 
we first expand it in terms of y-matrices and obtain a scalar, a pseudo scalar, a 
vector, an axial-vector and a tensor part. We then analyze the Lorentz structure of 
each part by expressing it in terms of possible “basic Lorentz covariants” and scalar 
functions. These scalar functions are known as TMD PDFs. There are totally 32 
such TMD PDFs. Among them, 8 contribute at leading twist and they all have clear 
probability interpretations such as the number density, the helicity distribution, the 
transversity, the Sivers function, the Boer-Mulders function etc; 16 contribute at 
twist 3 and the other 8 contribute at twist 4. We emphasize that they are all scalar 
functions of x and k±, i.e., depending on x and k^. If we integrate over (Pk±, terms 
where the basic Lorentz covariants are space odd vanish. At the leading twist, only 
3 of 8 survive, i.e. the number density, the helicity distribution and the transversity. 

Higher twist TMD PDFs are also defined via quark-j-gluon-quark correlators. 
Many of them are however not independent since they are related to those defined 
via the qu ark- quark correlator through the QCD equation of motion, ft is interesting 
to see thalPSl, although not generally proved, all the twist 3 TMD PDFs that are 
defined via quark-gluon-quark correlator and are involved in SlDlS are replaced by 
those defined via quark-quark correlator. 

1 also want to emphasize that fragmentation is just conjugate to parton distri¬ 
bution. We have one to one correspondence between TMD PDFs and TMD FFs. 


4. Accessing the TMDs in High Energy Reactions 

The TMDs can be studied in semi-inclusive high energy reactions such as SIDIS 
e~+N —>■ e~+h+X, semi-inclusive Drell-Yan h+h —>• /+-|-/“-|-A, and semi-inclusive 
hadron production in e+e“-annihilation e’*' -|- e“ hi + h 2 + X. With SIDIS, we 
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study TMD PDFs and TMD FFs, while with Drell-Yan and e“''e“ annihilation, we 
study TMD PDFs and TMD FFs separately. We now follow the same steps as those 
for inclusive DIS and briefly summarize what we already have in constructing the 
corresponding theoretical framework. 

(I) The general forms of hadronic tensors: For all three classes of processes, the 
general forms of hadronic tensors have been studied and obtained. For SIDIS, it has 
been discussed in Refs. [17ll^ and it has been shown that one need 18 independent 
structure functions for spinless h. For Drell-Yan, a comprehensive study was made 
in Ref. [21] and the number of independent structure functions is 48 for hadrons 
with spin 1/2. For e+e”-annihilation, the study was presented in Ref. [22] and one 
needs 72 for spin-1/2 hi and / 12 . 

(II) LO in pQCD and leading twist parton model results: These are the simplest 
parton model results and can be obtained easily. E.g., for SlDlS, the result can be 
obtained from those given e.g. in Ref. [20] by neglecting all the power suppressed 
contributions. 1 emphasize that the result obtained this way is a complete parton 
model result at LO in pQCD and leading twist. It can be used to extract the TMDs 
at this order. Any attempt to go beyond LO in pQCD or to consider higher twists 
needs to go beyond this expression. 

(III) LO in pQCD, leading and higher twist parton model results: For the semi- 
inclusive processes where only on e h adron is involved, either in the initial or the 
final state, it has been showiJMTj 

that the collinear expansion can be applied. 
Such processes include: semi-inclusive DIS e“ -I- —>■ e“ -I- q{jet) A, and e+e“- 
annihilation e"*" -I- e“ h q{jet) + X. By applying the collinear expansion, we 
have constructed the theoretical frameworks for these processes with which leading 
as well as higher twist contributions can be calculated in a systematical way to LO 
in pQCD. The complete results up to twist-3 have been obtained in Refs. P51l?7] .For 
unpolarized e~ + N ^ e~ + q{jet) -\- A, the results up to twist 4 have also been 
obtainecPSI. These results can be used as the basis for measuring these TMDs via 
the corres pon ding process at the LO in pQCD. I call in particular the attention to 
the results^ for + e~ ^ h + qijet) + X for ft, with different spins. Those for 
spin-1 hadrons involve tensor polarization that is mu ch less explored till now. See 
also talks by Y.K. Song and S.Y. Wei for more detailJ^SlMI 

However, for the above-mentioned three kinds of semi-inclusive processes, there 
are always two hadrons involved. Collinear expansion has not been proved how to 
apply for such processes. It is unclear how one can calculate leading and higher 
twist contributions in a systematical w ay. Never theless, twist 3 calculations that 
have been carried out for these processeJ^^EUSSj^ practically in the following steps: 
(i) draw Feynman diagrams with multiple gluon scattering to the order of one gluon 
exchange, (ii) insert the gauge link in the correlator wherever needed to make it 
gauge invariant, (iii) carry out calculations to the order 1/Q. Although not proved, 
it is interesting to see that the results obtained this way reduce exactly to those 
obtained in the corresponding simplified cases where collinear expansion is applied 
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if we take the corresponding fragmentation functions as (5-functions. 

(IV) TMD factorization and evolution: TMD factorization theorem has been 
established at the leading twist for semi-inclusive processei^S-lMl TMD evolution 
theory is also developing very fast^SHHO! See the overview talk by Daniel BoeiEH. 


5. TMD Parameterizations 


Experiments have been carried out for all three kinds of semi-inclusive reactions. 
The results are summarized in particular in the plenary talks in this conference by 
Marcin Stolarski, Renee Fatemi, and Armine Rostomyan. Here, 1 will just try to 
sort out the TMD parameterizations that we already have. 

The first part concerns what people called “the first phase parameterizations”, 
i.e. TMD parameterizations without QCD evolutions. 

(1) Transverse momentum dependence: This is usually taken aJ^^^^^a Gaussian 
in a factorized form independent of the longitudinal variable z or x. The width has 
been fitted, the form and flavor dependence etc have been tested. Roughly speaking, 
this is a quite satisfactory fit. However, it has also been pointed out, e.g. in [5S] for 
the TMD FF, that the Gaussian form seems to depend on the flavor and even on 
z, which means that it is only a zeroth order approxi mat ion. 

(2) Sivers function: It is usually parameterized^^^H^ in the form of the number 
density fq{x,k±) multiplied by a cc-dependent factor Afq{x) and a /c_L-dependent 
factor h{k±), and Nq{x) ^ x“(l — xY while h{kx) is taken as a Gaussian. There 
exist already different sets such as the Bochum, the Torino and the Vogelsang-Yuan 
fits. One thing seems to be clear that the Sivers function is nonzero for proton and 
it has different signs for u- and d-quark. 

(3) Transversity and Gollins function: A simultan eous e xtraction of them from 
SIDIS data have been carried out by the Torino groufJ^SlMl A similar form as that 
for the Sivers function has been taken and it has been obtained that also the Collins 
function is nonzero and has different signs e.g. for u tt’’' or d —>■ 7r+. 

(4) Boer-Mulders function: Clear signature for non-zero Boer-Mulders function 
has been obtained from SIDIS data on (cos 20) asymmetrji^Sl-Eni The form was 
taken again similar to the Sivers function. However, I would like to point out t hat 
the (cos 20) asymmetry receives twist-4 contributions due to the Cahn effedP^. 
A proper treatment of such twist-4 effect involves twist-4 TMDs as discussed in 
Ref. [24]. Because of the multiple gluon scattering shown in Fig. 1, the twist-4 
effects could be very much different from that given in [14] . 

Attempts to parameterize other TMDs such as pretzelocity have also been 
madi^. Although there is no enough data to give high accuracy constraints, the 
qualitative features obtained are also interesting. 

The second part concerns the TMD evolution. As mentioned earlier, this is a 
topic that develops very fast recently. A partial list of recent dedicated publications 
is Refs. [401150] . QGD evolution equations have been constructed. The numerical 
results obtained show clearly that TMD evolution is quite significant and it is im- 
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portant to use the comprehensive TMD evolution rather than a separate evolution 
of the transverse and longitudinal dependences respectively. See the overview talk 
by Daniel BoeiEB at this conference. 

At last, I want to mention the first version of TMD PDFlib has already 

created^. 

6. Summary and Outlook 

In summary, I just want to emphasize that three dimensional imaging of the nucleon 
is a hot and fast developing topic in last years. Many progresses have been made and 
many questions are open. Especially in view of the running and planned facilities 
such as the electron-ion colliders, we expect even rapid development in next years. 

I apologize for many aspects that I could not cover such as TMDs and Wigner 
function, model calculations, nuclear dependences, and hyperon polarization. The 
readers are referred to many interesting talks at this conference. 

I thank X.N. Wang, Y.K. Song, S.Y. Wei, K.B. Chen, J.H. Gao and many other 
people for collaboration and help in preparing this talk. My sincere thanks also go 
to John Collins for communications. This work was supported in part by NSFC 
(Nos.11035003 and 11375104), and the Major State Basic Research Development 
Program in China (No. 2014CB845406). 
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